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ABSTRACT 

We use a set of high-resolution simulations of scale-free Einstein-de Sitter cosmolo- 
gies to investigate the logarithmic slope of the phase-space density profile Q{r) — 
p{r)/a^{r) of dark matter (DM) haloes. The initial conditions for the simulations are 
determined by a power law power spectrum of the form P{k) cx fc". We compute the 
Q{r) profiles using the radial, tangential and full velocity dispersion, and the velocity 
anisotropy parameter, f3{r). We express Q{r) as a single power-law Q{r) oc r" and 
derive a median slope a in each simulation and for each definition of Q. Our main 
findings are: 1. The various Q{r) profiles follow a power law to a good approximation. 
2. The slopes depend on the concentration parameter c of the DM haloes, where for 
c > 10 the slopes steepen with rising concentration and for c < 10 the trend flattens 
and even turns around. 3. The asympt otic va l ue of (3 as r — > i?vir increases with the 
value of c. 4. In accordance with Zait et al. (20071 a^ad becomes more negative as 



the asymptotic value of (3 at the virial radius increases. 5. This introduces a weak 
dependence of the Q{r) slopes on the slope of the power spectrum. 

Key words: cosmology: theory — dark matter — methods: numerical 
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1 INTRODUCTION 

For the density profile p(r) of dark matter (DM) haloes a 
variety of empirical fitting formulae exist. The one applied 



mostly is the so-called NFW (Navarro et al. 1996 19971 



profile, which gives a fair description for haloes found in a 
wide range of A'^-body simulations. Recent studies though 
are hinting to use a different functional form for the den- 
sity profile that differs from the usual NFW form mainly 
in the central regions (see, e.g. Navarro et al. 2004 Mer 



ritt et al. 20061, but common to all suggested profiles is 



the fact that they are not described by a simple power-law 
but rather require a smoothly changing logarithmic slope. It 
comes therefore as a surprise that the coarse-grained phase- 
space densitjj^ profile, calculated from the density- and the 
velocity dispersion profile a{r) as 



Qir) 



p{r) 



(1) 



very closely follows a single power-law Q{r) oc with a ~ 
-1.9. This was first noted bylTaylor & Navarrol (120011) and 



has since been confirmed by many other works (e.g. Boylan- 



[Kolchin fc Ma|2004| |Rasia et al.|2"004l|Ascasibar et al.|20'0^ 



^ See, e.g., |Dehnen &: McLaughlm| ( [2005| for clarifications about 
the applicability of the term 'phase-space density', as the term is 
commonly used, however, we will continue to refer to Q{r) as the 
phase-space density profile. 



Austin et al.|[2005l [Dehnen fc McLaughlin] [20051 |Hoffman| 
et al.|2007||Ascasibar fc Gottlober|2008| ). 

Different velocity dispersions have been used with equa- 
tion [ij generally the radial velocity dispersion (Jiad('") is 
used, but also the total velocity dispersion atot{r) has been 
employed. Both lead to separate definitions for Q that nev- 
ertheless can be well fitted with a single power-law, albeit 



marginal differences in the slopes (Dehnen & McLaughlin 



2005||Faltenbacher et al.|2007||Ascasibar fc Gottlober|2008[ ), 



(20071 



see also the discussion in Zait et al 

The power law phase-space density profile has been 
found mostly in simulations of structure formation in the 
standard model of cosmology, namely the fiat ACDM model. 
It is not known how general its behaviour in other cosmolo- 
gies is and the dynamical origin of the power-law is not yet 
understood. However, the structure of (spherical) haloes in 
equilibrium obeys the (spherical) Jeans equation which re- 
lates the density profile, the phase-space density profile and 
the velocity anisotropy parameter /3 (to be defined below). 
As mentioned above, the universality of the density profile is 
well-established across a large variety of cosmological mod- 
els. The description of the phase-space density profiles as a 
single power-law, however, is not that well checked. 

Hence we aim here at studying the (phase-space) struc- 
ture of DM haloes in different cosmogonies, in particu- 
lar scale-free Einstein-de Sitter cosmologies with varying 
spectral index n, to explicitly probe universality of the 
phase-space profile. We use the numerical simulations of 
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Table 1. Simulation details. 



Name 


Type 


B 


A/as 






N 




512-0.50 


n = -0.50 


1 


0.0358 


1.0 


0.0 


5123 


186 


512-1.50 


n = -1.50 


1 


0.0215 


1.0 


0.0 


512^ 


106 


512-2.25 


n = -2.25 


1 


0.0109 


1.0 


0.0 


512^ 


61 


512-2.50 


n = -2.50 


1 


0.0076 


1.0 


0.0 


512^ 


34 


512-2.75 


n = -2.75 


1 


0.0046 


1.0 


0.0 


5123 


14 


256-B50 


ACDM 


50 


0.9 


0.3 


0.7 


256^ 


19 



Summary of the simulations. The first column simply assigns a 
unique label used throughout the paper while the second column 
specifies the model. B refers to the side length of the simulation 
box. A (erg) is the normalisation of the input power spectrum 
used to generate the initial conditions. Qm and Oa describe the 
background cosmology. N gives the number of particles in the 
model and the number of haloes used for the analysis. 



Knollmann et al. ( 2008 1 who studied the structure of haloes 
in pure DM flat cosmologies with a primordial power spec- 
trum of the form P{k) oc k", finding the equilibrium haloes 
to be well described by a universal density profile. Scale- 
free models provide a clean probe of the possible depen- 
dence of the structure of DM halos on the primordial initial 
conditions. This provides insight to the concordance ACDM 
model in the sense that the power spectrum for this cosmol- 
ogy exhibits a varying slope depending on the scale. 

We shall further investigate the behaviour of the phase- 
space density profile with different components of the total 
velocity dispersion crtot, namely the radial and tangential 
dispersions. Given the close connection between the Q{r) 



20071, 



power law and the velocity anisotropy (Zait et al 
expressed by the /3 parameter, the (3 profile of the DM haloes 
is to be explicitly studied. The paper is organized as follows. 
The numerical experiments are briefly discussed in §2. The 
fitting of the phase-space density profile is described in §3 
and the results are presented in §4. The paper concludes 
with a summary and discussion in §5. 



2 NUMERICAL EXPERIMENTS 

We base our analysis on the simulations described in |Knoll^ 



mann et al. (2008j). Those simulations are set up according 
to a power spectrum of the initial density perturbations of 
the form 



P{k) oc fc" 



(2) 



Five different models have been simulated with n = —0.50, 
— 1.50, —2.25, —2.50 and —2.75. The simulations were per- 



formed with the parallel A''-body code Gadget2 (Springel 
2005) following the evolution of 512'^ dark matter particles. 



The final stage of each run is chosen in such a way that M* 
- the typical collapsing mass - reaches the same value for all 



choices of n (for a more detailed description see Knollmann 
et al.|2008| |. 



We then employ the MPI version of the AM IGA Halo 
Findef ] (AHF, successor of MHF introduced by [Gill et al. 
( 2004[ )) to identify haloes and to compute their radial and 



integral properties. The radial properties are calculated in 
logarithmic bins in r and only those bins which are con- 
verged according to the criteria presented in [Power et al.| 
( 2003) are retained for the analysis. We further restrict our 
sample of haloes to those having a mass Al in the range 
[0.75M», 1.5M*] where a typical M* halo is resolved with 
4.2 X 10* particles. 

Also, we require the haloes to be relaxed, for which we 
have used the displacement Ar of the centre of mass of all 
material inside the virial radius with respect to the potential 
centre normalized to the virial radius 



(3) 



For inclusion in the sample, we require Ar < 0.05, employing 
a more conservative value than proposed by other recent 
studies (e.g. |Neto et al.|2007 l. 

For comparison we also use the results of a "standard" 
ACDM simulation analysed at z = whose particulars are 
described in detail elsewhere (Power et al. in prep.). Using 
that simulation we define a halo sample labelled 256-B50 
including all haloes resolved with more than 3.15 x 10* par- 
ticles (corresponding to a halo of mass 0.75M* in the scale- 
free simulations) and being relaxed according to the same 
"off-centre" criterion alluded to above. A summary of the 
simulation parameters and the halo samples is given in ta- 
ble [T] 



3 PHASE-SPACE DENSITY PROFILES 

The evaluation of the phase-space density as defined by 
equation [T] requires the (total) velocity dispersion atot (r) 
at the radius r, which is calculated as 



2 



(r) 



K 



Vtot,i 



(Vtot)l 



(4) 



where the sum extends over the A'^sh particles contained in a 
shell with radius r. With Vtot.i and (vtot) we refer to the ve- 
locity of the ith particle and the average velocity within the 
shell, respectively. For the calculation of the radial phase- 
space density Q^adir) we project the 3D velocity Vtot,i onto 
the radial direction and use only the radial component with 
equation |4j respectively. With these two velocity dispersions 
we can obtain the tangential velocity dispersion atan by us- 
ing the relation 



(5) 



for the calculation of the tangential phase-space density pro- 
file Qtan('-)- 



3.1 Normalization 

For each halo we normalize the radius r by r2 and the Q 
profile by Q2 = Q(r2), where r2 is the radius at which the 
density profile has a slope of — 2[^We find this "scaling ra- 
dius" in a non-parametric way by locating the maximum of 
r^p(r). Having identified r2 for each halo and representing 



^ AMIGA is freely available for download at |http : //mm . alp7| 
Ide/People/AKnebe/AMfGA/l 



3 In the case of an NFW profile r2 corresponds to the scale radius 
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its Q-profile with a cubic spline, Q2 is then readily available 
as the evaluation of the spline interpolation at r2. 

Additionally, we define for each halo a concentration 
parameter c as 



r-2 



(6) 



where r^ir is the virial radius defined via 3M(< rvir) = 
47rAvirP6T'vir with pi, being the cosmic background density 
and Avir=f78 (340) for the scale-free (ACDM) model. The 
concentration c has been found to be a distinguishing quan- 
tity between haloes from different scale-free models, with 
haloes forming in a shallower cosmology being more concen- 
trated (cf. ,Knollmann et al.|2008J . 



3.2 Fitting 

The normalized Q-profiles of each halo are fitted by a single 
power-law 



q(r;a) = Q2 



r2 



by minimizing the quantity 



Q{ri) - q{ri-a) 



■Q{rr) 



(7) 



(8) 



Note that due to the choice of our normalization we have 
only the slope a as a free parameter. 



4 RESULTS 

In this section, we will describe our results, starting in sub- 
section |4.f I with the power-law nature of the phase-space 
density profiles. We will then, in subsection |4.2[ focus on 
the individual slopes probing for the generality of the slopes 
and on what quantities they depend. We will then shortly 
investigate the relation between the three different Q pro- 
files in subsection |4.3[ before finally focussing on the velocity 
anisotropy profiles in subsection |4.4[ 



4.1 Phase-space density profiles 

In figure[T]the different Q-profiles (i.e. radial, tangential and 
total) are shown in the normalization described above along- 
side the best-fitting power-laws for the respective simulation 
(i.e. the power law as given by the median of all q's for a 
given model and Q, respectively). The median slopes used 
to draw the power laws are given in table [2] which we will 
describe later in greater detail. 

Furthermore, to avoid crowding, the profiles for the dif- 
ferent simulations are offset to one another by factors given 
in the key. For all profiles the power law structure is quite 
convincing - irrespective of the simulation. We also investi- 
gated the residuals of fits (though not shown here), finding 
the average deviation per degree of freedom to be on the 
order of 40 per cent, further strengthening the single power 
law fit be a reasonable choice. In the next subsection, we 
will use the slopes derived for each halo. 



4.2 The universality of the slope 

It is known that haloes formed in runs with different spectral 
index n favour different concentration^ but are otherwise 
compatible with a universal density slope (cf. |Knollmann| 



et al.j2008) 



We therefore investigate the slopes a of indi- 
vidual haloes as a function of their concentration c. This is 
shown in figure[2] where the left (middle, right) panel corre- 
sponds to the radial (angular, total) Q profile. Additionally, 
the simulations the haloes originate from are symbolised by 
different colours (see the caption of figure[2]for a description, 
we will use those symbols and colours for the remainder of 
the paper). The information contained in figure |2] is hence 
two-fold: a) the connection of a to the modeQ is visualised 
by the colour and b) the dependence on the concentration 
is encoded on the abscissa. 

To investigate those two effects more closely - and to 
understand which is the more fundamental dependence - 
we will probe the relations separately, looking for the one 
yielding the smaller scatter. To this extent, we bin all haloes 
from the scale-free models in c, regardless of their model 
identity; the medians and the upper and lower quartiles in 
the bins are shown in the left panel of figure [3] We then 
use the median a in each model to construct the right panel 
of figure |3j showing the trend for a — n. The curves for 
the different definitions of Q are offset to one another to 
avoid crowding. The median values and the upper and lower 
quartiles are also given in tables [2] and [3] respectively. 

We first note that, ignoring the model identity and 
combing all haloes, a clear, yet noisy, trend can be observed: 
For haloes with a concentration c > 10 the Q profiles steepen 
with rising c, whereas a fiattening of the a — c curve occurs 
for concentrations c < 10. This is reflected in the a — n 
relation, for flat n models, we flnd steeper Q proflles, with 
the n = —2.75 model, however, showing a steeper median a 
than the n — —2.50 and n = —2.25 models. We also note 
that neither the a — c nor the a — n is clearly preferred as 
the fundamental relation, both have a comparable scatter in 
a, with maybe the former slightly favoured. 

In order to check whether our results are an artifact of 
the fltting procedure we performed three additional test: 
a) we fltted the Q-profiles confining the radial range to 
0.1 < r/r2 < 10 (i.e. the fitting has been performed over 
a fixed dynamical range for all the haloes of all the differ- 
ent models considered here), b) we subdivided the sample 
of each model into a low and high concentration bin and 
calculated the median a in those two bins for each model, 
and c) we calculated the reduced for each individual halo 
fit and inspected its distribution. All these checks indicate 
that the results seen and presented in figures [2] and |3] are 
stable. 

For test a) we recover a similar distribution as shown 
in figure [2] The second test traces the global trend shown in 
the left panel of figure[3] with the n = —0.50 and n — —1.50 
models showing a steeper Q-profile for the high concentra- 
tion subsample, whereas the n — —2.25 shows a similar a 



■* This is already apparent in figure [T] the Q-profiles are plot- 
ted out to the virial radius and hence the final point is a direct 
measure of concentration for our choice of units. 
^ This means effectively a dependence on the spectral index n of 
the power spectrum of initial density perturbations. 
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10'^ 10° 10^ 10"'' 10° 10' lO"' 10° io' 

r/rg r/rj r/rj 



Figure 1. Phase-space density profiles Q(r) for our sample of haloes (radial, tangential, and total velocity dispersions are used in the 
left, middle, and right panel). The profiles of the different simulations are shifted on the y-axis for clarity, the factor by which they are 
shifted is given in the key. For each simulation we also show the best-fitting r" power-law, as given by the median value of a within the 
simulation; cf. table [2] for the actual values. 



Table 2. Median slopes of the phase-space density profiles de- 
pending on the model. 



Model 



512-0.50 
512-1.50 
512-2.25 
512-2.50 
512-2.75 
256-B50 



-2.02: 
-1.951 
-1.891 
-1.891 
-1.9i: 
-1.94" 



2.08 
1.88 
2.02 
1.82 
1.98 
1.82 
1.93 
1.82 
1.94 
1.83 
1.98 



-2.oo: 

-1.94; 
-1.92; 

-1.87: 

-1.92: 
-1.81 



-1.91 
-2.06 
-1.83 
-2.02 
-1.86 
-2.02 
-1.79 
-1.93 
-1.84 
-2.02 
-1.75 
-1.84 



-1.99: 
-1.92: 
-i.9o: 
-1.87: 
-1.93: 

-1.85 



-1.93 
-2.04 
-1.86 
-2.00 
-1.86 
-1.99 
-1.83 
-1.92 
-1.83 
-1.99 
-1.82 
-1.87 



Table 3. Median slopes of the phase-space density profiles of the 
scale-free haloes depending on the concentration. 



"rad 



5 056.20 
'-'■""-'3.69 


-1.921 


1.82 
1.99 


-1.931 


1.86 

2.07 


-1.92 


q cic;10.83 
"•"-''-'8.34 


-1.921 


1.86 
2.02 


-1.901 


1.78 
1.98 


-1.91 




-1.951 


1.89 
2.00 


-1.911 


1.83 
1.97 


-1.90 


Qo 01 39.61 
00. ZI27. 81 


-1.981 


1.94 
2.05 


-1.991 


1.91 
2.03 


-1.98 


C7 cr:64.01 
■^'■"^''51.00 


-2.031 


2.00 
2.08 


-2.041 


1.99 
2.09 


-2.03 



-1.85 
-2.02 
-1.84 
-1.93 
-1.87 
-1.95 
-1.93 
-2.02 
-1.99 
-2.07 



in both bins and the steeper models show a steeper a for 
the low concentration bin, e.g. the opposite of what is found 
in the flatter models. It should be noted though, that due 
to the low numbers of haloes available in the steep models 
(n < —2.25, cf. table [1]), this might just be a statistical ef- 
fect. Lastly, the probability distributions of xLdi Xtan, and 
Xtot (not shown here though) confirm that all profiles are 
fitted equally well with the assumption of a single power- 
law, with the relative deviations per degree of freedom of 
the order of 40 per cent. 



4.3 The relation between the different Q profiles 

We will now investigate the relation of the different Q pro- 



McLaughlin] |2005[ [Ascasibar fc Gottlober] |2008[ ) we find 
Qrad to be steeper than Qtot in our ACDM sample, i.e. 
(f^rad/«tot)ACDM ~ 1.05. The scale-frcc simulations show 
at best a very marginal trend of this kind: The slopes of 
the radial and total Q-profiles for a given model tend to be 
much less different, i.e. (Q!r/Qtot)scaio-free ~ 1.01. However, 
the slope for a given Q-profile appears to change across mod- 
els, namely to drop with increasing spectral index n (shown 
in the right panel of [3|, with the notable exception of the 
512-2.75 model which favours steeper Q profiles than the 
512-2.50 model. As we have shown above, this is also seen 
as a a — c dependence, namely the slope of the Q profiles 
first rises with decreasing concentration c until c « 10, from 
where on the trend fiattens and even turns around; this is 
shown in the left panel of fi gure [3] 



Ascasibar & Gottlober ( 2008 1 were the first to point out 



that not only Qrad and Qtot are single power-laws, but also 
Qtan- They further showed that its slope is shallower than 
the slope Qtot basing their findings on a suite of simulations 
of the concordance ACDM cosmology. We confirm this re- 
lation for the haloes in the (concordance) 256-B50 sample; 
we further find a single power-law to be a fair description 
of the Qtan profile in the scale-free samples. However, in the 
scale-free sample atan tends to be slightly steeper than Oftot, 
but still shallower than Orad. 

4.4 The velocity anisotropy of scale-free haloes 

Given the universal density profile the difference between 
the power law slopes of the different Q(r) profiles implies 
that the anisotropy (3 of the velocity dispersion, given as 



1 



^2 
"tan 



(9) 



files. In agreement with previous studies (e.g. Dehnen & 



is not constant. We investigate this in figure [4] showing in 
the left panel the mean j3{r) profiles of the different models. 
This shows an agreement of all profiles for r < r2 but a 
dependence on the model in the outer part of the haloes. 
In the middle and right panel we now focus on the value 
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-1.6 
-1.7 
-1.8 
-1.9 
B -2.0 
-2.1 
-2.2 
-2.3 
-2.4 



-'rad 





^tot 




10 



10 



10 



Figure 2. In this figure we show the derived slopes a of each halo versus its concentration c = c^-^jr^. Each panel corresponds to using 
a different velocity dispersions to calculate Q(r') (see equationjlj, the results for Qrad('') {Qtan(^'), QtotirY} are shown in the left (middle, 
right) column. The haloes are colour coded according to the model they correspond to: red n = —0.50 (open boxes), green n = —1.50 
(open circles), blue n = —2.25 (open upright triangles), magenta n = —2.50 (open downright triangles), cyan n = —2.75 (open diamonds) 
and for comparison haloes from a ACDM simulation are shown in orange (open pentagons). 



-1.4 



-1.5 



-1.6 



-1.7 



-1.9 



-2.0 



-2.1 



radH-0.4 



tanH-0.2 



4-4 



tot 



4 4" 



I ' ' ' ' I 




t 



10" 



10' 



-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 
n 



Figure 3. In this figure we investigate the relation of the slope a of the phase-space density profiles to the concentration c (left panel) 
and the spectral index n of the model (right panel). For clarity, the points corresponding to the tangential (dark green filled circles) and 
radial (dark red filled squares) Q-profile are offset to the ones corresponding to the total (dark blue filled triangles) Q-profile by 0.2 and 
0.4, respectively. In the left panel we have binned the a — c distribution shown in figure[2]in c regardless of n, whereas in the right panel 
we combine the slopes of all haloes within one model regardless of the concentration (see also tables [2] and [sjl . We show the median and 
the lower and upper quartiles of each bin. 



of the f5 profile at the virial radius and check for how it 
correlates with the slope of the radial Q profile Qrad and the 
concentration c, respectively. 

We see a trend for /3(rvir) to become larger for a steeper 
radial Q profile, which also translates into a trend with 
c, namely, the larger the concentration the larger becomes 
/3(f'vir); this is shown in the right-most panel of figure |4] 



Recalling equation [9] we see that a larger /3(rvir) means 
that the total velocity dispersion becomes more radial. The 
/3(rvir) — Orad trend we observe agrees with the findings of 



Zait et al. ( 2007 1 who calculated the f3 profile given the den- 
sity profile and a power law Qrad profile. They showed that 
a more negative a^ad implies a larger (3 for r > r2, which 
means that the total velocity dispersion is more radial. 
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Figure 4. In the left panel, we show the mean (3 profile in each simulation, using the same colour coding as described in figure |2] Note 
that the last point corresponds to value of (3 at the virial radius. The middle panel gives the value of the anisotropy profile at the virial 
radius, /3(rvir), as a function of the slope of the radial Q-profile, Qrad- In the right panel the distribution of /3{rvir) as a function of the 
concentration is shown. 



5 SUMMARY & DISCUSSION 

We have analysed five high-resolution simulations of scale- 
free cosmologies and investigated the logarithmic slope a 
of the radial, angular and total phase-space density pro- 
files Q{r) = p{r)/a^(r) when fitted by a single power-law 
Q{r) oc r" . Additionally, we used a halo sample drawn from 
a ACDM simulation as a control sample to verify our anal- 
ysis procedure and reproduced previous results. 

Focussing on a well-defined sample of haloes in each 
scale-free simulation, chosen to be in a comparable dynami- 
cal state, we derived the slopes aiad, cetan and atot for each 
halo and computed the median slopes in each simulation. We 
found that the slopes for the three definitions of the phase- 
space density profile Q vary between the models, namely Q 
becomes shallower with steeper spectral index n of the initial 
power spectrum P{k) oc fc", with the exception of the 512- 
2.75 model, which shows fiatter Q profiles than found in the 
512-2.50 model. We also find a trend for a to first become 
flatter with decreasing concentration; this trend flattens at 
c ~ 10 and even turns around for the tangential and total 
Q proflles. This is connected to previous flndings showing 
a correlation between the power spectrum index n and the 
concentration c. 

But how can the n - c relation affect the Q{r) pro- 
flles? The scaled density profiles appear to follow a univer- 
sal form, independent of the cosmological model (see e.g. 
Navarro et al.|1996 1997 2004 1 ; in particular this has been 
shown explicitly for the models studied here ( |Knollmann| 
et al.|2008 K The dependence on the cosmological model (i.e. 
the slope of the power spectrum of initial density perturba- 
tions) is introduced via the concentration parameter, and 
thereby by the dependence of the /3 profile on the model. 
The Jeans equation, which dictates the structure of spher- 
ical DM haloes in equilibrium, relates the density, phase 



which is precisely what we find and have shown in figure |4] 
It is hence interesting to note that, even though the den- 
sity profile follows a universal form, the /9-profile does not, 



space density and the velocity anisotropy profiles (e.g. Zait 
et al.|[2007 and references therein). Zait et al. showed that 
for a given density profile a more radial velocity dispersion 
at the outer parts of a halo implies a more negative arad. 



but rather depends on the concentration (Knollmann et al 
[20081). 



Summarizing past and present results the following con- 
clusions can be drawn about DM haloes in scale-free cos- 
mologies. The density profile shows a universal profile, upon 
an appropriate scaling, but the concentration parameter de- 



pends on the power spectrum index n (Knollmann et al, 
[2008). Here it has been found the phase-space density pro- 
file obeys a power law. The slope of its power law depends 
on whether the phase-space is defined in terms of the ra- 
dial, tangential or the full velocity dispersion. Furthermore, 
the slopes depends on the shape of the power spectrum. 
Also, the velocity anisotropy at i?vir depends strongly on 
the power spectrum index. 



The main problem addressed here is the universality of 
the phase-space density profile, and a mixed answer has been 
found. On the one hand it has been very robustly shown that 
indeed a power law profile provides a good fit in all models 
considered here. Yet, the power law slope a has been found 
to vary with power spectrum index n. Given our lack of 
understanding of the origin of the phase-space density power 
law it is interesting to trace the reason for the dependence 
of the slopes a on n. It has been shown here that at least 
part of that dependence is attributed to the concentration 
parameter variation with n. However, it is unclear whether 
there is a direct dependence of the slope a on n. This can 
be tested by comparing the slope and n relation for a sub- 
sample of haloes of different models with the same value 
of the concentration parameter. Regrettably, the sample of 
haloes used in the present work is too small to provide a 
clear answer. 
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